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We report a measurement of CP asymmetry parameters in the decay B''{B^) —* J/tp-K^ , which is 
governed by the 6 — » ccd transition. The analysis is based on a 140 fb~^ data sample accumulated 
at the T(4S') resonance by the Belle detector at the KEKB asymmetric-energy e^e~ collider. We 
fully reconstruct one neutral B meson in the J/ipTv^ final state. The accompanying B meson flavor 
is identified by its decay products. From the distribution of proper time intervals between the two B 
decays, we obtain the following CP-violating parameters: Sj/^^o = — 0.72 ± 0.42(stat) ±0.09(syst) 
and Aj/^^o = -0.01 ± 0.29(stat) ± 0.03(syst). 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw 



In the Standard Model (SM), the Kobayashi-Maskawa 
(KM) quark-mixing matrix ij has an irreducible complex 
phase that gives rise to CP violation in weak interac- 
tions. In particular, the SM predicts large CP-violating 
asymmetries in the time-dependent rates of P'^ and 
decays into a common CP eigenstate fcp jH- In the de- 
cay chain T(45) — > P"P° /cp/tag, where one of the 
B mesons decays at time tcp to a final state fcp and 
the other decays at time ttag to a final state /tag that 
distinguishes between P° and P°, the decay rate has a 
time dependenct probability given by Q 



P(At) = 



4tro 



1 + g 



sin(AmdAi) 



(1) 



where tbo is the P° lifetime, Anid is the mass difference 
between the two P° mass eigenstates. At = tcp ~ ^tag, 
and the 6-flavor q — +1 (—1) when the tagging B meson 
is a P" (P'^). The CP-violating parameters Sf^,p and 
Afcp are given by 
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IAp + 1' 
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where A is a complex parameter that depends on both 
the B'^B^ mixing and on the amplitudes for P" and P° 
decay to fcp- To a good approximation in the SM, |A| is 
equal to the absolute value of the ratio of the B'^ ^ fcp 
to B'^ ^ fcp decay amplitudes. 

CP violation in neutral P meson decays involving the 
h — > CCS transition has been established through mea- 
surements of the CP- violating parameter sin 2(j)i by the 



Belle 13] and BaBar _5j collaborations. The SM predicts 
Sfcp = — i^/sin20i, where ^/ — +1(— 1) corresponds to 
CP-even (-odd) final states; and -A/^p — (or equiv- 
alently |A| = 1) for tree diagrams in both b — > ccs and 
6 — > ccd. In contrast with the h ccs case, however, tree 
and penguin amplitudes all contribute to 6 — > ccd to the 
same order in the sine of the Cabibbo angle. Therefore, 
if penguin or other contributions are substantial, a preci- 
sion measurement of the time-dependent CP asymmetry 
in 6 ^ ccd may reveal values for 5j/^^o and Aj /^,t^o that 
differ from the values for h ccs. Measurements of CP 
asymmetries in 6 — > ccd transition-induced P decays such 
as B'^ —> J/f/jTT^ thus play an important role in probing 
these one-loop diagrams. 

A study of CP asymmetry in P*^ J/i/jtt^ decays has 
been reported by the BaBar collaboration based on 
81 fb~^. In this paper we report a measurement of time- 
dependent CP-violating parameters in P° J/ipir'^ de- 
cays using a data sample of 140 fb~^ (corresponding to 
15.2x10^ PP pairs) collected at the T(4S') resonance 
with the Belle detector at the KEKB asymmetric- 
energy e+e" (3.5 GeV on 8.0 GeV) coUider (ij. The 
T(4iS') is produced with a Lorentz boost of /37 — 0.425 
nearly along the z-axis defined as anti-parallel to the 
positron beam. Since the B^ and P'^ mesons are nearly 
at rest in the T(4S') center-of-mass system (cms), At 
can be determined from Az, the displacement in z be- 
tween the fcp and /tag decay vertices: At ~ Az/Pjc = 

[ZCP - Ztag)//37C 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a three-layer silicon vertex de- 
tector (SVD), a 50-layer central drift chamber (CDC), an 
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array of aerogel threshold Cherenkov counters (ACC), 
a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a su- 
perconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the 
coil is instrumented to detect mesons and to iden- 
tify muons (KLM). The selection of hadronic events is 
described elsewhere . 

J/i/) mesons are reconstructed via their decays into 
oppositely-charged lepton pairs (e+e" or iS^ Both 
lepton tracks must be positively identified [9j. In 
the e+e" mode, ECL clusters that have no associated 
charged tracks and are within 50 mrad of the track's 
initial momentum vector are included in the calcula- 
tion of the invariant mass (Mee(-^)), in order to in- 
clude photons radiated from electrons/positrons. The 
invariant masses of e^e^{'^) and /i^/i^ combinations 



are required to fall within the ranges 
(Mee(^) - Mji^) < 36 MeV/c^ and 



-150 MeV/c^ < 



sured properties of these tracks, two parameters, q and r, 
are assigned to each event. The first, g, has the discrete 
value -1-1 (—1) when the tag-side B meson is more likely 
to be a {B'^). The parameter r is an event-by-event 
MC-determined flavor-tagging quality factor that ranges 
from r = for no flavor discrimination to r = 1 for an un- 
ambiguous flavor assignment. It is used only to sort data 
into six intervals of r, according to the estimated flavor 
purity. The wrong-tag probabilities, wi and the differ- 
ence in wrong-tag probabilities between 5° and S", Awi 
{I = 1,6) are fixed using a data sample of self-tagged 5° 
decay modes. The wrong tag fractions for each r interval 
that are used in the final fit are given elsewhere • 

The decay vertices of B mesons are reconstructed using 
tracks that have associated SVD hits. Each vertex posi- 
tion is required to be consistent with the interaction point 
profile smeared by the average transverse B meson decay 
length 12]. The vertex position for the B^ J/tpir'^ 
decay is reconstructed using lepton tracks from the J/ip- 



respectively. Here M, 



J/i: 



(M^^ - Mj/^) < 36 MeV/c2, ^ „ 
denotes the world average of the J/ 'ijj mass jl^ • 

Photon candidates are selected from clusters of up to 
5x5 crystals in the ECL. Each photon candidate is re- 
quired to have no associated charged track, and a cluster 
shape that is consistent with an electromagnetic shower. 
To select tt*^ 77 decay candidates for the — > J/ tpir'^ 
mode, the energy of each photon is required to ex- 
ceed 50 MeV in the ECL barrel (32° < < 129°) or 
100 MeV in the forward (17° < 9 < 32°) or backward 
(129° < 6 < 150°) endcaps, where 9 denotes the polar 
angle with respect to the z-axis. Neutral pion candidates 
are formed from photon pairs that have an invariant mass 
in the range 0.118 GeV/c^ to 0.150 GeV/c^. 

J/ijj and tt'^ candidates are combined to form B can- 
didates. The B candidate selection is carried out using 
two observables in the T(4S') cms: the beam-energy con- 
strained mass Mbc = V-E-beam ~ iJ^Pi)'^ the energy 
difference AE = Y,Ei- i^beam, where -Bbeam = is 
the cms beam energy, and pi and Ei are the cms three- 
momenta and energies of the B meson decay products, re- 
spectively. In this calculation, pi and Ei are obtained af- 
ter refitting with vertex and mass constraints for the J /if) 
di-lepton decays and a mass constraint for the tt*^ — > 77 
decays in order to improve the AE' and A/bc resolutions. 
The B meson signal region is defined as 5.27 GeV/c^ < 
Mbc < 5-29 GeV/c2 and -0.10 GeV < A£; < 0.05 GeV. 
The lower bound of AE is chosen to accommodate the 
negative AE tail due to shower leakage associated with 
7r°, and avoid the background at AE ^ —0.2 GeV due 
to B^ J/i/) Ks{Ks 7r'^7r°) events. The number of 
reconstructed 5° J/ipir'^ candidates is 103. 

Charged leptons, kaons, pions, and A baryons that are 
not associated with the reconstructed B'^ J/-07r° can- 
didate are used to identify the 6-flavor of the accompany- 
ing B meson, which decays into /tag- Based on the mea- 



-60 MeV/c < The /tag vertex is determined from all well- reconstructed 



tracks except those used for B° J/ipir^ reconstruc- 
tion and tracks that form a Kg candidate. After fla- 
vor tagging and vertex reconstruction, 91 of the 103 
B'^ J/ijjTT^ candidates remain. The AE and Mbc dis- 
tributions for the candidate events are shown in Fig. ^ 




-0.2 -0.1 0.1 
AE (GeV) 



FIG. 1; (a) Ai5 distribution for events in the Mbc signal re- 
gion and (b) A/bc distribution for events in the /S.E signal re- 
gion, for J/ip-K^ candidates. The superimposed curves 
show fitted contributions from signal (dot-dashed line, plot- 
ted only in Ai5 distribution) , B J/^l)X background (dotted 
line) , combinatorial background (dashed line) and the sum of 
all the contributions (solid line). 

To assign an event-by-event signal probability for 
use in the maximum-likelihood fit to extract the CP- 
violating parameters, we determine event distribution 
functions in the Ai?-Mbc plane for both signal and back- 
ground. The signal distribution is modeled with a two- 
dimensional function which is Gaussian in Mbc and a 
Crystal Ball line shape [3| in AE. Since the number 
of candidate events is not sufficient to obtain precisely 
the shape parameters of these functions, those are de- 
termined from MC simulation and held fixed in the fit. 
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while the overall signal yield is allowed to float. The 
background distribution is studied using a large sample 
of MC events along with events outside the signal region. 
We split the backgrounds into two categories: B decays 
containing a J/t/j {B — > J/ipX) and combinatorial back- 
ground to which random combinations of particles in BB 
decays and continuum events contribute. According to 
a MC study, the B J/tpX background forms a small 
peak in the Mbc projection. Therefore we parametrize 
this contribution with a sum of a Gaussian and a phase- 
space-like background function (ARGUS function) ll^l in 
Mbc and an exponential function for AE. The amount 
of i? — > J/ tpX background is estimated by the MC study 
and fixed in the fit For the combinatorial back- 

ground, we use a linear function for AE and an AR- 
GUS function for Mbc. From the fit, the purity of the 
B'^ J/tpTT^ candidates is estimated to be (86±11)%. 

We determine Sjj^T^o and Aj j^T^a by performing an un- 
binned maximum-likelihood fit to the observed At distri- 
bution. The probability density function (PDF) expected 
for the signal distribution is given by 



^sig(At) 

g-IAtl/rgO 



4t 



so 



\-qAwi^q{\~2wi) 



Sji^^o sin(AmrfAt) -f Aji^,^o cos(AmrfAt) 



(3) 



to account for the effect of incorrect flavor assignment. 
The distribution is convolved with the proper-time in- 
terval resolution function i?sig(Ai), which takes into ac- 
count the finite vertex resolution. i?sig(At) is formed 
by convolving four components: the detector resolutions 
for zcp and ztag, the shift in the Ztag vertex position 
due to secondary tracks originating from charmed par- 
ticle decays, and the kinematic approximation that the 
B mesons are at rest in the cms j^. A small compo- 
nent of broad outliers in the Az distribution, caused by 
misreconstruction, is represented by a Gaussian function 
Pol (At). We use the same resolution function parameters 
as used in the sin20i measurement llj. We determine 
the following likelihood value for each event: 

Pl{Ati\Sj/^T,0,Aj/^T,o) 



(1-/01) j_ 



/,ig7'sig(At')i?sig(At, - At') 



bkg '-bkg mKLt i^u-At') 



+-'bkg ' 1 ^ 

+ (l-/sig-/b/g^'')nkg'(^0 



xi?™f (At,-At') 



a!(At') + /oiPoi(At,), (4) 



where /oi is the outlier fraction, /sig and /b/^"^ are the 
signal and B J/il}X background probabilities calcu- 

be, Viit''{At) and 



V^'^^iAt) are the PDFs for B J/-0X and combinato- 
rial background events, respectively. V^^^^{At) is mod- 
eled as a sum of exponential and prompt components, 
and is convolved with the corresponding resolution func- 
tion -Rbkg*' which is modeled by a sum of two Gaussians. 
The parameters in V{^'^{At) and are determined 

by a fit to the At distribution of a background-enhanced 
control sample, i.e. events not in the Ai?-Mbc signal re- 
gion. V^l^'^ {At) is a lifetime distribution determined by 
MC simulation. The contributions of CP = — 1 and +1 
components in the B — s- J /ipX background in the signal 
region are found to be almost the same and the CP vio- 
lation effect is largely cancelled. We model -Rb^g using 
the signal resolution function because vertex reconstruc- 
tion for both signal and B — > J/ipX background events 
is based on lepton tracks from J/^. We fix r^o and Am^ 
at their world-average values [ig . The only free param- 
eters in the final fit are Sj/^^^o and .4j/^^o, which are 
determined by maximizing the likelihood function 



(5) 



where the product runs over all events. 

A fit to the candidate events results in the CP- 
violating parameters, 

Sj/^^o = -0.72±0.42(stat)±0.09(syst) 
Aj/^^o = -0.01±0.29(stat)±0.03(syst), (6) 

where the sources of systematic errors are described be- 
low. Figure |2 (a) and (b) show the At distributions for 
W J/V'tt" (g = +1) and P° ^ J/i/'7r° (g = -1) event 
samples respectively. Figure |21 (c) shows the raw asym- 
metry in each At bin without background subtraction, 
which is defined by 



Nq= + l + Ng=-l ' 



(7) 



lated as a function of AE and M 



where Nq=+i{Nq=^i) is the number of observed candi- 
dates with q = +1(— 1). The curve shows the result of 
the unbinncd maximum-likelihood fit to the At distribu- 
tions. 

The major systematic errors come from uncertainties 
in the vertex reconstruction (±0.06 for Sj/^^^a and ±0.02 
for Aj/^jro), the wrong tag fraction (±0.03 and ±0.01), 
the potential CP- violating effect in i? ^ J/i/jX back- 
ground (±0.03 and ±0.01), the signal probability (±0.02 
and ±0.02), and possible fit bias (±0.03 and ±0.01). The 
contributions from other systematic error sources, un- 
certainties in the background At distribution, resolution 
function, and physics parameters {tbo, Anid), are found 
to be much smaller. The quadratic sum of all the con- 
tributions amounts to ±0.09 for Sj/^-^a and ±0.03 for 

Aj/tPt^o. 

In summary, we have performed a measurement of 
the CP- violating parameters in the P° J/tpn'^ decay. 
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FIG. 2: The At distributions for (a) B° J/tPtv° (q = +1) 
and (b) B° J/ip-K^ {q = -1) candidates, (c) The asym- 
metry, A, in each At bin. The curve shows the result of the 
unbinned maximum- likelihood fit, and dashed curves show 
the background distributions. 



The resultant values are S 



J I'ip-n 



0.09(syst) and Aj/^^o = -0.01 ± 0.29(stat) ± 0.03(syst). 
These values are consistent with those obtained for 
J jil) Ks and other decays governed by 5 — > ccs transition, 
as expected from the SM if the tree diagram dominates. 
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